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Abstract. In this paper, we study quasi-stationarity for a large class of Kolmogorov dif- 
fusions. The main novelty here is that we allow the drift to go to — oo at the origin, and 
the diffusion to have an entrance boundary at +00. These diffusions arise as images, by a 
deterministic map, of generalized Feller diffusions, which themselves are obtained as limits 
of rescaled birth-death processes. Generalized Feller diffusions take nonnegative values and 
are absorbed at zero in finite time with probability 1. An important example is the logistic 
Feller diffusion. 

We give sufficient conditions on the drift near and near +00 for the existence of quasi- 
stationary distributions, as well as rate of convergence in the Yaglom limit and existence 
of the Q-process. We also show that under these conditions, there is exactly one quasi- 
stationary distribution, and that this distribution attracts all initial distributions under the 
conditional evolution, if and only if +00 is an entrance boundary. In particular this gives 
a sufficient condition for the uniqueness of quasi-stationary distributions. In the proofs 
spectral theory plays an important role on L 2 of the reference measure for the killed process. 
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1. Introduction 

The main motivation of this work is the existence, uniqueness and domain of attraction of 
quasi-stationary distributions for some diffusion models arising from population dynamics. 
After a change of variable, the problem is stated in the framework of Kolmogorov diffusion 
processes with a drift behaving like —l/2x near the origin. Hence, we shall study quasi- 
stationarity for the larger class of one-dimensional Kolmogorov diffusions (drifted Brownian 
motions), with drift possibly exploding at the origin. 

Consider a one-dimensional drifted Brownian motion on (0, 00) 

dX t = dB t - q(X t ) dt , X = x > (1.1) 

where q is defined and C 1 on (0, 00) and (Bf]t > 0) is a standard one-dimensional Brownian 
motion. In particular q is allowed to explode at the origin. A pathwise unique solution of 
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(jl.ip exists up to the explosion time r. We denote T y the first time the process hits y G (0, oo) 
(see |14j chapter VI section 3) before the explosion 

T y = inf{0 < t < t : X t = y} 

We denote by = lim T n and To = lim Ti/ n . Since q is regular in (0, oo) then r = ToAT^. 

n— +oo n^oo ' 

The law of the process starting from Xq with distribution v will be denoted by ¥ u . A quasi- 
stationary distribution (in short q.s.d.) for X is a probability measure v supported on (0, oo) 
satisfying for all t > 

P„(Xi G 4 | T > t) = u(A), V Borel set A C (0, oo). (1.2) 

By definition a g.s.d. is a fixed point of the conditional evolution. The Yaglom limit tt is 
defined as the limit in distribution 

tt(.) = lim F x (X t G • | T > t), 

t—*oo 

provided this limit exists and is independent of the initial condition x. The Yaglom limit is 
a q.s.d. (see Lemma [772]) . 

We will also study the existence of the so-called Q-process which is obtained as the law of 
the process X conditioned to be never extinct, and it is defined as follows. For any s > 
and any Borel set B C C([0, s]) consider 

Q X (X G B) = lim ¥ X (X £B\T >t). 

t— >oo 

When it exists, this limit procedure defines the law of a diffusion that never reaches 0. 

The reason for studying such diffusion processes with a possibly exploding drift at the origin 
comes from our interest in the following generalized Feller diffusion processes 



dZ t = ^JWtdB t + h{Z t )dt, Z = z>0, (1.3) 
where h is a nice function satisfying h(0) = 0. 

Notice that z = is an absorbing state for Z. This means that if Zq = then Zt = for all 
t is the unique solution of (jl.3p (see fT4"]). 

If we define X t = 2^Z t /j then 

dX t = dB t -^- (---h (?^X) dt , X = x = l^fzff > 0, (1.4) 



X t \2 7 V 4 

so that X is a drifted Brownian motion as in (|l.lj) where g(a;) behaves like l/2x near the 
origin. The process Z is obtained after rescaling some sequences of birth-death processes 
arising from population dynamics. 

Of particular interest is the case h(z) = rz — cz 2 (logistic case), for which we obtain 

1 rx c^fx 3 

A complete description of these models is performed in the final section of the paper, where 
their biological meaning is also discussed. Of course quasi-stationary distributions for Z 
and X are related by an immediate change of variables, so that the results on X can be 
immediately translated to results on Z. 
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The study of quasi-stationarity is a long standing problem (see [29] for a regularly updated 
extensive bibliography and [101 113|, [33] for the Markov chain case). For Kolmogorov diffusions, 
the theory started with Mandl's paper [26J in 1961, and was then developed by many authors 
(see in particular [271 135]), All these works assume Mandl's conditions, which are not 
satisfied in the situation described above, since in particular the drift is not bounded near 

0. It is worth noticing that the behavior of q at infinity also may violate Mandl's conditions, 
since in the logistic case for instance, 

QW f e Q{y) dydz < oo, 
Ji Ji 

where Q(y) := 2 J® q(x) dx. 

This unusual situation prevents us from using earlier results on q.s.d.'s of solutions of Kol- 
mogorov equations. Hence we are led to develop new techniques to cope with this situation. 
In Section [2] we start with the study of a general Kolmogorov diffusion process on the half 
line and introduce the hypothesis (HI) that ensures to reach in finite time with probability 

1. Then we introduce the measure fj,, not necessarily finite, defined as 

fx(dy) := e~ Q M dy, 

which is the speed measure of X. We describe the Girsanov transform and show how to 
use it in order to obtain L 2 (^) estimates for the heat kernel (Theorem I2.3|) . The key is the 
following: starting from any x > 0, the law of the process at time t is absolutely continuous 
with respect to /_i with a density belonging to L 2 (/j) (and explicit bounds). In the present 
paper we work in L p (/j) spaces rather than \?{dx), since it greatly simplifies the presentation 
of the spectral theory. 

This spectral theory is done in Section [3l where we introduce the hypothesis (H2): 
lim q 2 (x) — q'(x) = oo , C := — inf q 2 (x) — q'(x) < oo. 

This hypothesis ensures the discreteness of the spectrum (Theorem I3.2p . The ground state 
r/i (eigenfunction associated to the bottom of the spectrum) can be chosen nonnegative, even 
positive as we will see, and furnishes the natural candidate T)\d\i for a q.s.d.. The only thing 
to check is that rji £ L 1 (/-t) which is not immediate since \i is possibly unbounded. 
Section [J] gives some sharper properties of the eigenfunctions defined in the previous section, 
using in particular properties of the Dirichlet heat kernel. More specifically, we introduce 
two independent hypotheses, (H3) and (H4), either of which ensures that the eigenfunctions 
belong to L x (/j) (Propositions 14.31 and 14. 4p . Hypothesis (H3) is 

fi(dy) < oo, 



q 2 (y)-q'(y) + C + 2 



and (H4) is 



e Q ^ dy < oo , ye 

Q(y)/2 dy 

< oo. 

J Jo 



Section [5] contains the proofs of the existence of the Yaglom limit (Theorem I5.2f) as well 
as the exponential decay to equilibrium (Proposition 15. 5h . under hypotheses (HI) and (H2), 
together with either (H3) or (H4). Section[6]contains the results on the Q-process (Corollaries 
EJ] and E2]). 
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In Section [7] we introduce condition (H5) which is equivalent to the existence of an entrance 
law at +00, that is, the repelling force at infinity imposes to the process starting from 
infinity to reach any finite interval in finite time. The process is then said to "come down 
from infinity" . We show that the process comes down from infinity if and only if there exists 
a unique q.s.d. which attracts any initial law under the conditional evolution (Theorem 17.3ft . 
In particular this theorem gives sufficient conditions for uniqueness of q.s.d.'s. In the context 
of birth and death chains the equivalence between uniqueness of a q.s.d. and "come down 
from infinity" has been proved in [8j Theorem 3.2]. 

The final section contains the description of the underlying biological models, as well as the 
application of the whole theory developed in the previous sections to these models (Theorem 
E2}. 

In the following statement (which is basically Theorem 18. 2ft , we record the main results of 
this paper in terms of the generalized Feller diffusion solution of (11.30 . We say that h satisfies 
the condition (HH) if 

/ ■ \ r Kx) xh'(x) 

m lim — — = —00, (11) lim = 0. 

x^oo yjx z->oo h(X) z 

Theorem 1.1. If h satisfies (HH), then for all initial laws with bounded support, the law of 
Zt conditioned on {Zt 7^ 0} converges exponentially fast to a probability measure v , called the 
Yaglom limit. 

The process Z conditioned to be never extinct is well-defined and is called the Q-process. 
The Q-process converges in distribution, to its unique invariant probability measure. This 
probability measure is absolutely continuous w.r.t. v with a nondecreasing Radon-Nikodym 
derivative. 

If in addition, the following integrability condition is satisfied 

dx 

1 ~h(x) 

then Z comes down from infinity and the convergence of the conditional one- dimensional 
distributions holds for all initial laws, so that the Yaglom limit v is the unique quasi- stationary 
distribution. 

2. One dimensional diffusion processes on the positive half line 

Associated to q we consider the functions 

A(s) = j* e Q{y) dy and k{x) = f* e Q{y) e~ Q{z) dz^j dy , (2.1) 

where we recall that Q(y) = Jf 2q{u)du. Notice that A is the scale function for X. 
For most of the results in this paper we shall assume sure absorption at zero, that is 
Hypothesis (HI): for all x > , F x (t = T < T^) = 1 . (2.2) 
It is well known (see e.g. [2] chapter VI Theorem 3.2) that (|2.2p holds if and only if 

A(oo) = 00 and k(0 + ) < cc . (2.3) 

We notice that (HI) can be written as ¥ x ( lim A tAr = 0) = 1. 

t— >oo 
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Example 2.1. The main cases that we are interested in are the following ones. 

(1) When X is defined by (jl.4p associated to the generalized Feller diffusion Z. It is 
direct to show that Q(x) behaves like \og(x) near hence k(0 + ) < oo. The logistic 
Feller diffusion corresponds to h(z) = rz — cz 2 for some constants c and r. It is easily 
seen that (|2,3|) is satisfied in this case provided c > or c = and r < 0. 

(2) When the drift is bounded near 0, in which case k(0 + ) < oo. 



We shall now discuss some properties of the law of X up to To. The first result is a Girsanov 
type result. 

Proposition 2.2. Assume (HI). For any Borel bounded function F : C([0, t], (0, oo)) — ► R 
it holds 



E x [F{X) l t<To ] 



K" 



F (v) WoH exp ( ]- Q(x) - ]- Q(u t ) 



(q 2 ~ q')(v s )ds 



where W™ x denotes the expectation w.r.t. the Wiener measure starting from x , andE x denotes 
the expectation with respect to the law of X starting also from x. 

Proof. It is enough to show the result for F nonnegative and bounded. Let x > and 
consider e 6 (0, 1) such that e < x < 1/e. Also we define r e = T £ A Ti/ £ . Choose some ip £ 
which is a nonnegative C°° function with compact support, included in ]e/2,2/s[ such that 
ip e ( u ) = life<^<l/e. The law of the diffusion (II. ip coincides up to t £ with the law of a 
similar diffusion process X £ obtained by replacing q with the cutoff q £ = qvp £ . For the latter 
we may apply the Novikov criterion ensuring that the law of X £ is given via the Girsanov 
formula. Hence 



E x [F(X) I 



t<T E 



E 



w x 



F (u) It< Te (w) ex P (J -qe(Us)duJs ~\ J G 



(q £ ) 2 {u s )ds 



E 



w x 



F(uj) l t <r e {w) exp 



E 



w x 



(uj s )dw s 
1 



I t<r e (^) exp Q Q(x) - ^ Q(ut) ~ \ J 



(u) s )ds 



(q 2 - q')(u s )ds 



The last equality is obtained integrating by parts the stochastic integral. But !f< Te is non- 
decreasing in e and converges almost surely to lt<T both for ¥ x (thanks to (HI)) and W^. 
It remains to use Lebesgue monotone convergence theorem to finish the proof. □ 



The next theorem is inspired by the calculation in Theorem 3.2.7 of [32]. It will be useful to 
introduce the following measure defined on (0, oo) 

H(dy) := e~ Q{y) dy . 

Note that /i is not necessarily finite. 



(2.4) 



Theorem 2.3. Assume (HI). For all x > and all t > there exists some density r(t,x, .) 
that satisfies 



E x [f(X t ) l t<To ] 



f(y)r(t,x,y)n(dy) 



for all bounded Borel f . 



6 P. CATTIAUX, P. COLLET, A. LAMBERT, S. MARTINEZ, S. MELEARD, AND J. SAN MARTIN 

If in addition there exists some C > such that q 2 (y) — q'{y) > — C for all y > 0, then for 
all t > and all x > 0, 



r 2 (t,x,y)fi{dy) < (l/27rt)* e Ct e Q{ ^ . 



Proof. Define 



G(u) = l t <T (u) exp Q Q(u ) - ~ Q(u t ) ~ \ J O? 2 ~ q'){^s)ds^j 



Denote by 



-v(t,x,y) 



(2nt) 2 exp 



2t 



the density at time t of the Brownian motion starting from x. According to Proposition 
we have 

E x [f(X t )l t<To ] = E W *[/(^)E W *[GM] 

/(y)E w ^[GK = y] e~ v ^ dy 

f(y) E w * [G\«h = y] e -*C^.w)+OC») , 

because E W:c [G|tJt = y] = if y < 0. In other words, the law of Xt restricted to non-extinction 
has a density with respect to [i given by 

r(t, x, y) = E w * [G\uh = v] e - v{t ' x ' y)+Qiy) . 

Hence, 



r 2 (t,x,y) fi(dy) 



E* 



< E* 



E w *[G\uj t = y] e^M+flW) e -Q(v)+v(tp,v) e -«(*^.v) d y 



e -i)(t,x,o; t )+Q(w t ) ]gW x [C 2 !^] 



o(t,x,u t ) p - Jo {q 2 -q'){u s )ds 



where we have used Cauchy-Schwarz's inequality. Since e v ^^ x '^ < (l/27ri)2 the proof is 
completed. □ 

Remark 2.4. It is interesting to discuss a little bit the conditions we have introduced. 

(1) Since q is assumed to be regular, the condition q 2 — q' bounded from below has to be 
checked only near infinity or near 0. 

(2) Consider the behavior near infinity. Let us show that if liminf J/ _^ OC) (g 2 (y) — q'(y)) = 
— oo then limsup J/ ^ 00 (g 2 (y) — q'(y)) > — oo i.e. the drift q is strongly oscillating. 
Indeed, assume that q 2 {y) — q'(y) — * — oo as y — » oo. It follows that q'(y) — ► oo, hence 
q(y) — > oo. For y large enough we may thus write q(y) = e u ^ for some u going to 
infinity at infinity. So e 2u<y ' y \l — v! \y)e~ u ^) — * — oo implying that u'e~ u > 1 near 
infinity. Thus if g = e~ u we have g' < —1 i.e. g(y) — > —oo as y — > oo which is 
impossible since g is nonnegative. 
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(3) If X is given by (jl.4p we have 




Hence, since h is of class C 1 and h(0) = 0, q 2 {y) — q'{y) behaves near like ^ so that 

q 2 — q' is bounded from below near (see Appendix for further conditions fulfilled by 
h to get the same result near oo). 



3. L 2 and spectral theory of the diffusion process 

Theorem 12.31 shows that for a large family of initial laws, the distribution of Xt before 
extinction has a density belonging to L 2 (/i). The measure fi is natural since the kernel of the 
killed process is symmetric in L 2 (/i), which allows us to use spectral theory. 
Let Co°((0, oo)) be the vector space of infinitely differentiable functions on (0, oo) with com- 
pact support. We denote 

poo 

(f,9)ix = / f(u)g(u)fi(du) . 
Jo 

Consider the symmetric form 

£(/,<?) = </',<?% , = C °°((0,oo)). (3.1) 

This form is Markovian and closable. The proof of the latter assertion is similar to the one of 
Theorem 2.1.4 in [11] just replacing the real line by the positive half line. Its smallest closed 
extension, again denoted by £, is thus a Dirichlet form which is actually regular and local. 
According to the theory of Dirichlet forms (see [11] or [12]) we thus know that 

• there exists a non-positive self adjoint operator L on L 2 (^x) with domain D(L) D 
C(5°((0, oo) such that for all / and g in Cq°((0, oo)) the following holds (see [TT] 
Theorem 1.3.1) 

/•oo 

£(f, g) = -2 / f(u) Lg{u) ^{du) = - 2 (/, Lg) „ . (3.2) 
Jo 

We point out that for g £ C^°((0, oo)), 

Lg = \g"-qg>. 

• L is the generator of a strongly continuous symmetric semigroup of contractions on 
L 2 (//) denoted by {Pt)t>o- This semigroup is (sub)-Markovian, i.e. < Ptf < 1 fi 
a.e. if < / < 1 (see [33] Theorem 1.4.1). 

• There exists a unique ^t-symmetric Hunt process with continuous sample paths (i.e. 
a diffusion process) up to its explosion time r whose Dirichlet form is £ (see |11] 
Theorem 6.2.2) 

The last assertion implies that, for /i quasi all x > (that is, except for a set of zero capacity, 
see [H] for details), one can find a probability measure Q x on C(M + , (0,oo)) such that for 
aU/€C °°((0,oo)), 

/(%r) - f(x) - / Lf(uj s )ds 
Jo 
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is a local martingale with quadratic variation Jq At \f'\ 2 (u> s )ds. Due to our hypothesis q G 
C 1 (0, oo) we know that this martingale problem admits a unique solution (see for example 
|15] page 444). On the other hand, using Ito's formula we know that under P x , the law of 
(Xt/\ T ) is also a solution to this martingale problem. 

The conclusion is that the semigroup Pt and the semigroup induced by the strong Markov 
process (X t /\ T ) coincide on the set of smooth and compactly supported functions. Therefore, 
for all / G L 2 (m) we have that 

P t f(x) = E[f(X t )l t<T ]. 



Let (E\ : A > 0) be the spectral family of —L. We can restrict ourselves to the case A > 
because — L is nonnegative. Then Vi > 0, /, g G h 2 (/j,), 

/roo 
P t fgdn = J e- xt d(E x f,g)^. (3.3) 

We notice that if absorption is sure, that is (HI) holds, this semigroup coincides with the 
semigroup of X killed at 0, that is P t f(x) = E[/(X?) l t <T ]- 
Note that for / G L 2 (//) and all closed interval K C (0, oo), 

J{P t f?dn = J(P t (fl K + fl K c)) 2 dfi 

< 2 J (P t (fl K )) 2 d f l + 2 j (P t (fl K r)) 2 d^i 

< 2 J (P t (fl K )) 2 d f i + 2 J (flKofdfi. 

We may choose K large enough in order that the second term in the latter sum is bounded 
by e. Similarly we may approximate fix in L 2 (^) by fix for some continuous and bounded 
/, up to e (uniformly in t). Now, thanks to (HI) we know that Pt(/l^)(x) goes to as t 
goes to infinity for any x. Since 



J {P t (fl K )) 2 dfi = fP 2 t(fl K ) d^i . 



we may apply Lebesgue bounded convergence theorem and conclude that J (Pt(flK)) 2 dfi — > 
as t — > oo. Hence, we have shown that, 

V / G 7L 2 (n) J (Ptffdfi -> as t -» oo . (3.4) 
Now we shall introduce the main assumption on q for the spectral aspect of the study. 



Hypothesis (H2) C = — inf q 2 (y) — q'{y) < oo and lim q 2 (y) — q'{y) = +oo . (3.5) 

3/S(0,oo) y^co 

Proposition 3.1. Under (H2), \q(x)\ tends to infinity as x — > oo, and q~(x) or q + (x) tend 
to as x [ 0. If in addition (HI) holds then q(x) — ► oo, as x — > oo. 
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Proof. Since q 2 — q' tends to oo as x — ► oo, q does not change sign for large x. If q is 
bounded near infinity we arrive to a contradiction because q' tends to — oo and therefore q 
tends to — oo as well. So q is unbounded. If limfnf \q(x)\ = a < oo then we can construct 

x^oo 

a sequence x n — > oo of local maxima, or local minima of q whose value |g(a? n )| < a + 1, but 
then q 2 {x n ) — q'{x n ) stays bounded, which is a contradiction. 

Now we prove that q~{x) or q + (x) tend to as x [ 0. In fact, assume there exist an e > 
and a sequence (x n ) with < x n { such that q(x2n) = —e, q(x2n+i) = £• Then we can 
construct another sequence z n j such that |g(zn)| < e an d l'(z n ) — ► oo, contradicting (H2). 
Finally, assume (HI) holds. If g(x) < —1 for all x > x$ we arrive to a contradiction. Indeed, 
for all t 

^o+i(7o >t)> ^x Q +i{T X0 >t)> Px +i{Tx = oo) . 
The assumption q(x) < — 1 implies that Xt > Bt + t while t < T xo , and therefore 

e -2(xo+l) _ g 2a;o 

F X0+1 (T >t)> F xo+1 (B t + t hits oo before x ) = e _ 2oo _ e2a;o = 1 - e" 2 

where we have used that (exp(— 2(Bt + t))) is a martingale. This contradicts (HI) and we 
have q(x) — > oo as x — > oo. □ 



We may now state the following result. 

Theorem 3.2. If (H2) is satisfied, —L has a purely discrete spectrum < Ai < A 2 < 
Furthermore each Aj (i E f$) is associated to a unique (up to a multiplicative constant) 
eigenfunction rji of class C 2 ((0,oo)) ; which also satisfies the ODE 

1 

2*?/ ~ 1Vi = -\r)i- (3-6) 

The sequence (?7i)i>i is an orthonormal basis of L 2 (//) ; rji can be chosen to be strictly positive 

in (0, oo). 

ForgGlL 2 ^), 

Ptg = J2 e ~ Xit ( r >i>9)nm in h 2 (ji), 

then for f,g£ 

lim e Ali ( 5 ,P t /) M = (vi,f)n(m,9)n- 

t — >oo 

//, in addition, (HI) holds, then X\ > 0. 

Proof. For / G h 2 (dx), define P t (/) = e~ Q / 2 P t (f e Q l 2 ), which exists in L 2 ((ix) since / e Q l 2 G 
L 2 (^). (A)t>o is then a strongly continuous semigroup in h 2 (dx), whose generator L coincides 
on Cg°((0,oo)) with \£j - \ (q 2 - q') since C 2 (0,oo) C D(L), and e Q l 2 G C 2 (0,oo). The 
spectral theory of such a Schrodinger operator on the line (or the half line) is well known, but 
here the potential v = (q 2 — q')/2 does not necessarily belong to L^ c near as it is generally 
assumed. We shall use [3] chapter 2. 

First we follow the proof of Theorem 3.1 in [3] . Since we have assumed that v is bounded from 
below by — C/2, we may consider H = L — (C/2 + 1), i.e. replace vhyv + C/2 + l = w> 1, 
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hence translate the spectrum. Since for / G Co°(0,oo) 

/■OO /'OO /'OO 

-(Hf, /):=-/ Hf(u)f(u)du = {\f'{u)\ 2 /2 + w(u)f 2 {u))du> f(u)du, (3.7) 
jo io Jo 

H has a bounded inverse operator. Hence the spectrum of H (and then the one of L) will be 
discrete as soon as H^ 1 is a compact operator, i.e. as soon as M = {/ G D{H) ; —{Hf, f) < 
1} is relatively compact. This is shown in [3j when w is locally bounded, in particular 
bounded near 0. If w goes to infinity at 0, the situation is even better since our set M is 
included into the corresponding one with w ~ 1 near the origin, which is relatively compact 
thanks to the asymptotic behavior of v. The conclusion of Theorem 3.1 in [3] is thus still 
true in our situation, i.e. the spectrum is discrete. 

The discussion in Section 2.3 of [3], pp. 59-69, is only concerned with the asymptotic behavior 
(near infinity) of the solutions of /" — 2wf = 0. Nevertheless, the results there applies to 
our case. All eigenvalues of L are thus simple (Proposition 3.3 in [3]), and of course the 
corresponding set of normalized eigenfunctions (ipk)k>i is an orthonormal basis of h 2 (dx). 

The system (e^/ 2 ipk)k>i is thus an orthonormal basis of L 2 (//), each rjk = e^/ 2 ?/^ being an 
eigenfunction of L. We can choose them to be C 2 ((0, oo)) and they satisfy (|3.6j) . 
For every t > 0, and for every g, f G L 2 (/^) we have 

oo 

Y, e ~ Xkt ^g)^ k j)^ = {g,P t f)^. 
k=l 

In addition if g and / are nonnegative we get 

0< lim e Al *( 5 ,PJ) / , = (7 7l) /) M (r/ 1 , g ) / ,, 

t — >oo 

since Ai < A2 < . . . and the sum Y^kLi li^k, d) ^iVk, f)fi\ is finite. It follows that (t/i,/)^ 
and (r]i,g)^ have the same sign. Changing r/i into —r/i if necessary, we may assume that 
(Viif)fi ^ f° r an Y nonnegative /, hence 771 > 0. Since Ptrji(x) = e~ Xlt rji(x) and 771 is 
continuous and not trivial, we deduce that r]\{x) > for all x > 0. 

Since L is non-positive, Ai > 0. Now assume that (HI) holds. Using (|3.4p we get for g G L 2 (/i) 

= lim (P t g, P t g)^ = lim e- 2Xlt (g, r^) 2 , 

showing that Ai > 0. □ 

Moreover, we are able to obtain a pointwise representation of the density r. 

Proposition 3.3. Under (HI) and (H2) we have 

00 

r(t, x,y) = Y^ e- Xkt m(x)ri k (y), (3.8) 

k=l 

uniformly on compact sets o/(0,oo) x (0, 00) x (0,oo). 
Therefore on compact sets of (0, 00) x (0, 00) we get 

lim e Xlt r(t, x,y) =r ll {x)r}\ (y) . (3.9) 

t—fOO 
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Proof. Using Theorems 12.31 and 13.21 for every smooth function g compactly supported on 
(0, oo) we have 

n oo „ 

fc=i fc=i J J 

Then using the regularity of rj k and r we obtain, by letting g(y)dy tend to the Dirac measure 
in x, that 

n 

Y, e ' Xkt Vk{x) 2 <r{t,. 



. _ N ,2,2). 
fc=l 

oo 

Thus, the series ^ e~ Xkt rj k (x) 2 converges pointwise, which by Cauchy-Schwarz inequality 

k=l 

oo 

implies the pointwise absolute convergence of x,y) := ^ e~ Xktr q k (x)r] k (y) and the bound 

k=l 

for all n 



Afc *|%(2)%(y)| < y / r(t,x,x)y / r(t,y,y). 



k=l 



Using Harnack inequality (see for example [20]) we get 



y/r(t,x,x)^r(t,y,y) < C K r(t,x,y) 

for any x and y in the compact subset K of (0, oo). Using the dominated convergence theorem 
we obtain that for all Borel functions g, f with compact support in (0, oo) 

g{x)f(y)at,x,y)e-^-^dxdy= [ [ g(x)f(y)r(t,x,y)e~^-^dxdy. 



Therefore ((t,x,y) = r(t,x,y) dxdy-a.s., which proves the almost sure version of ()3.8f) . 
Since r\ k are smooth eigenfunctions we get the pointwise equality 

e- x ^r] k {x) 2 = e- Afc '/ 3 (r(t/3, x, •), Vk )^(r(t/3, x, •),ri k ) li 

= f J r(t/3,x,y)r(t/3,x,z)e- Xkt / 3 r ]k (y)ri k (z)e- Q(z " > - Q ^dydz, 

which together with the fact r(t/3, 2,») G L 2 (//) and Theorem 12.31 allow us to deduce 

ET=i e~ Afc V(2) 2 = / / r(t/3, x, y)r(t/3, x, z) ^ k=1 e^/^y^^e-Q^-Qiv) dydz 
= J J r(t/3, x, y)r(i/3, x, z)r(t/3, y, z)e~^^~^^ dydz = r(t, x, x). 

Dini's theorem then proves the uniform convergence in compacts of (0, oo) for the series 

oo 

^e~ Xkt r] k (x) 2 = r(t,x,x). 

k=l 

By the Cauchy-Schwarz inequality we have for any n 

/ oo \ 1/2 / oo \ 1/2 



J^e Xkt m(x)r] k (y) 



k=n 



i ' 'ik\y^ 2 

\k=n / \k=n 



This together with the dominated convergence theorem yields (|3.9p . □ 

In the previous theorem, notice that X^fc e_Afe * = / r (t,x,x)e~^^dx is the L 1 ^) norm of 
x i — ^ r(i, 2, x). This is finite if and only if Pt is a trace-class operator on L 2 (//). 
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4. Properties of the eigenfunctions 

In this section, we study some properties of the eigenfunctions rji, including their integrability 
with respect to fj,. 

Proposition 4.1. Assume that (HI) and (H2) are satisfied. Then r}ie~®dx < oo, F(x) = 
rj[(x)e~Q( x ^ is a nonnegative decreasing function and the following limits exist 

F(0 + ) =limr ? / 1 (x)e" Q(a:) G (0,ool, F(oo) = lim r][(x)e~ Q ^ G [0,oo). 

xlO x^oo 
oo , 

Moreover J r/i(x)e~^^dx = F ^°°^ . In particular 

o 1 

7/1 G if and only if F(0 + ) < oo. 

oo 

The function rji is increasing and J e~^ y 'dy < oo. 

l 

Remark 4.2. Note that g = r}\e~Q satisfies the adjoint equation ^g" + (qg)' = —X\g, and 
then F(x) = g'(x) + 2q(x)g(x) represents the flux at x. Then r/i G L 1 (/j) or equivalently 
g G L 1 ((ix) if and only if the flux at is finite. 

Proof. Since r/i satisfies rf{(x) — 2qr]' 1 (x) = — 2Air/i(x), we obtain for xq and x in (0, oo) 

rj[(x)e' Q ^ = rj[(x )e- Q ^ - 2Ai [* m{y)e~ Q{y) dy, (4.1) 

and F = n[e~® is decreasing. Integrating further gives 

m ( x ) = m (x ) + I ( r ] ' 1 (x )e~ Q ^ - 2Ai ( Vi(y)e~ Q{y) dyj e Q{z Uz. 

If for some zq > xo it holds that n' 1 (xo)e~Q( x °^ — 2Ai f*° r]i(y)e~Q( y 'dy < 0, then this inequal- 
ity holds for all z > zq since the quantity 

' 1 (x )e-^-2X 1 f Z m (y)e-^dy 



J XO 

is decreasing in z. This implies that for large x the function r\\ is negative, because e®^ 
tends to oo as z — > oo. This is a contradiction and we deduce that for all x > 

/•OO 

2Ai / m (y)e- Q ^dy <n[(x)e- Q ^. 

J X 

This implies that r\\ is increasing and, being nonnegative, it is bounded near 0. In particular, 

oo 

r/i(0 + ) exists. Also we deduce that F > and that J e~®^dy < oo. We can take the limit 

l 

as x — ► oo in (|4.ip to get 

F(oo) = lim 7/i(x)e~ Q ^ G [0,oo), 

ie— >oo 

and 77^ (xq)^^ ) = F(oo) + 2Ai ni(y)e~^^dy. From this equality the result follows. □ 
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In the next results we give some sufficient conditions, in terms of q, for the integrability of 
the eigenf unctions. A first useful condition is the following one 

Hypothesis (H3): / ^— — — — - e~ Q ^dy < oo, 

Jo i (y) - Q {y) + c + 2 

where as before C = — inf (q 2 (x) — q'(x)). 

Proposition 4.3. Assume that (HI), (H2) and (H3) are satisfied. Then rji belongs toL 1 (/x) 
for all i. 

Proof. Recall that ipi = e~®/ 2 r]i is an eigenfunction of the Schrodinger operator H introduced 
in the proof of Theorem 13.21 Replacing / by ipi in (|3.7|) thus yields 



roo poo 

/.'|2/„ \ /o i . \./.2/ 



m z (y)/2 + w(y)^(y))dy . 



poo pc 

(C/2 + 1 + A.) / $(y)dy= / 
Jo Jo 

Since the left hand side is finite, the right hand side is finite, in particular 

' poo 

w{y)n 2 {y)ii{dy) = \ w{y)ip 2 {y)dy < oo . 



As a consequence, using Cauchy-Schwarz inequality we get on one hand 



L 



i i 
, ..»-// i 



1 / rl \ a / r 1 1 \ 2 



Vi(y)\Kdy) < ( / w(y)n i (y)n(dy)) [ \ —^-^(dy)) < oo 



'0 \Jo 
thanks to (H3). On the other hand 



i 

2 * - / / \ - 



\rn(y)\n(dy) < (/ rj i (y)fj,(dy) 1 M*/)J < oo 

according to Proposition 14.11 We have thus proved that 77$ G L 1 (^). □ 

We now obtain sharper estimates using properties of the Dirichlet heat kernel. For this reason 
we introduce 



/oo ^1 
e -Q( x )dx < oo and J 



x e~^l 2 dx 



< oo. 



Proposition 4.4. Assume (H2) and (H4) hold. Then all eigenf unctions rjk belong to L 1 (//) ; 
and there is a constant K\ > such that for any x G (0, oo) and any k 

\Vk(x)\ <K X e Xk e Q{x)/2 . 

Moreover rji is strictly positive on M + , and there is a constant K 2 > such that for any 
x G (0, 1] and any k 

\Vk(x)\ <K 2 xe 2Xk e Q(x)/2 . 
Proof. In Section [3] we introduced the semigroup Pt associated with the Schrodinger equation 

Q 

and showed that r/k = e 2 tpk, where ipk is the unique eigenfunction related to the eigenvalue 
Afc for Pf. Using estimates on this semigroup, we will get some properties of ipk, and we will 
prove the proposition. 
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The semigroup p is given for / £ L 2 (R + , dx) by 
P t f(x) = E r 



f{u(t)) i t<To exp (-^ J o (i 2 ~ q'){u s )ds 



where E W:c denotes the expectation w.r.t. the Wiener measure starting from x. We first 
establish a basic estimate on its kernel pt(x,y). 

Lemma 4.5. Assume condition (H2) holds. There exists a constant K% > and a continuous 
increasing function B defined on [0,oo) satisfying lim 2 ^ 00 i?(z) = oo, such that for any 
x > 0, y > we have 

< pi{x,y) < e -(*~y?/± e -B(^{x,y}) ( 42 ) 

and 

p 1 (x,y) <K 3 p°(x,y) , (4.3) 
where p® is the Dirichlet heat kernel in R + given for M + by 

1 / -(x-y) 2 (x+y) 2 ' 



e 2t 



/2nt 

The proof of this lemma is postponed to the Appendix. 

It follows immediately from the previous lemma that the kernel pi(x,y) defines a bounded 
operator Pi from L 2 (K + , dx) to L°° dx) . As a byproduct, we get that all eigenfunctions 
iph of Pi are bounded, and more precisely 

\i/>k\ < K x e Afe . 

One also deduces from the previous lemma that the kernel defined for M > by 

vf{x,y) = l x< Mt y<M Pi{x,y) 

is a Hilbert-Schmidt operator in L 2 (M + , <ix) , in particular is a compact operator (see for 
example [7j pages 177, 267]). In addition, it follows at once again from Lemma 14.51 that if 
denotes the operator with kernel pf 1 , we have the following estimate, in the norm of 
operators acting on L 2 (M + , dx) , 

where C is a positive constant independent of M. Since ]imM- too B(M) = oo, the operator Pi 
is a limit in norm of compact operators in L 2 (R + , dx) and hence compact. Since pi(x, y) > 0, 
the operator Pi is positivity improving (that is if ^ f > then Pi/ > ) implying that 
the eigenvector tp\ is positive. 

We now claim that IV'fc^)! < K% x e 2Afc for < x < 1. We have from Lemma 14.51 and the 
explicit expression for p®(x,y) the existence of a constant K% such that 



e~ x 



oo / n /■oo 



k M*) <^ 3 y o P?(x,y)\Mv)\ dy < KsWfaW^ -e~ x 12 J ^ e~y I 2 S iifo{xy) dy . 
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We now estimate the integral in the right hand side. Using the convexity property of sinh 
we get sinh(xy) < xsinh(y) < |e y , for x £ [0, l],y > which yields 

roc „, roo 



roo poo 

/ e~ y2/2 sinh(xy) dy < - / e - y2/2 e y dy 
Jo 2 y 



proving the claim. Together with hypothesis (H4), this estimate implies that % belongs to 
^((0,1),^). 

Since 
we have 

/"OO /-oo 

r/ fc d/i = / ^fc(x) e~ Q ^ 2 dx 



which implies % £ L 1 ((l, oo), dfij using Cauchy-Schwarz's inequality. This finishes the proof 
of Proposition 14.41 □ 

Remark 4.6. Let us discuss some easy facts about the hypotheses introduced. 

(1) If q and q' extend continuously up to 0, hypotheses (H2), (H3) and (H4) reduce to 
their counterpart at infinity. 

(2) Consider q(x) = § + g{x) with g a C 1 function up to 0. In order that (|2.3p holds 
at the origin we need a > —5. Then fi(dx) = Q(x)x~ 2a dx with O bounded near the 
origin, while q 2 (x) — q'(x) ~ (a + a 2 )/x 2 . Hence for (H2) to hold, we need a > 0. 
Now we have the estimates 

' n{dx)^ [ £ @(x)x 2{1 - 

Jo 



L-a) 



/o (q 2 {x)-q'(x)+C + 2) 

and 



JO 

Therefore (H3) holds for a < § and (H4) holds (at 0) for a < 2. The conclusion is 
that a £ [0, f ). 

We recall that a = \ if X comes from a generalized Feller diffusion. 

(3) If q(x) > for x large, hypothesis (H2) implies the first part of hypothesis (H4). 
Indeed, take a > be such that for any x > a we have q(x) > and q 2 (x) — q'(x) > 1. 
Consider the function y = e~®l 2 which satisfies y' = —qy and y" = (q 2 — q')y. For 
b > a we get after integration by parts 

= / {(q 2 ~ q')y 2 - yy") dx = Uq 2 - q)y 2 + y' 2 ) dx - y(b)y'(b) + y[a)y'{a) . 
Using y' = —qy we obtain 

y 2 dx < J ({q 2 - q)y 2 + y' 2 ^j dx = q(a)y(a) 2 - q(b)y(b) 2 < q(a)y(a) 2 < 00 
and the result follows by letting b tend to infinity. 
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5. Quasi-stationary distribution and Yaglom limit 

Existence of the Yaglom limit and of q.s.d. for killed one-dimensional diffusion processes 
have already been proved by various authors, following the pioneering work by Mandl [26] 
(see e.g. [51 [271 [35] and references therein). One of the main assumptions in these papers is 
«(oo) = oo and 

e- Q ^ (J* e Q( - z) dz^j dy = +oo 

which is not necessarily satisfied in our case. Indeed, under mild conditions, the Laplace 
method yields that J? e@\ z ) dz behaves like e®^ /2q(y) when y tends to infinity, so the 
above equality will not hold if q grows too fast to infinity at infinity. Actually, we will 
be particularly interested in these cases (our forthcoming assumption (H5)), since they are 
exactly those when the diffusion "comes down from infinity", which ensures uniqueness of 
the q.s.d.. The second assumption in the aforementioned papers is that q is C 1 up to the 
origin which is not true in our case of interest. 

It is useful to introduce the following condition 

Definition 5.1. We say that Hypothesis (H) is satisfied if (HI) and (H2) hold, and more- 
over 7/i G L 1 ( / u) (which is the case for example under (H3) or (H4)). 

We now study the existence of q.s.d. and Yaglom limit in our framework. When r\\ £ L 1 (/j), 
a natural candidate for being a q.s.d. is the normalized measure 771/^/(771, which turns 
out to be the conditional limit distribution. 

Theorem 5.2. Assume that Hypothesis (H) holds. Then 

rfidfi 

is a quasi- stationary distribution, that is for every t > and any Borel subset A of (0, 00), 

F ul (X t eA\T >t) = vi(A). 
Also for any x > and any Borel subset A of (0, 00) , 

lime Al *P a .(r >t) = t/i(x)(t/i,1) m , (5.1) 



lim e Xlt F x (X t £ A, Tq > t) = v±{A) 7/1 (x) (7/1, 1) 



t—*oo 



l-l ■ 



This implies since 7/1 > on (0, 00) 

lim ¥ x (X t e A\T >t) = vx(A), 

and the probability measure v\ is the Yaglom limit distribution. Moreover, for any probability 
measure p with compact support in (0, 00) we have 



lim e Al * P P (T > t) = (771, 1)„ f m (x)p(dx) ; (5.2) 

lim e x ^F p (X t £A,T >t) = v\(A) (rji, 1)^ / m (x) p(dx) ; (5.3) 
lim FJX t £ A I T > t) = vi(A). (5.4) 
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Proof. Thanks to the symmetry of the semigroup, we have for all / in L 2 (/j), 

PtfTfidpL = J fPtmdfi = e~ Xlt J fvidfi. 

Since 771 £ L 1 (/i), this equality extends to all bounded /. In particular we may use it with 
/ = I(o, 00) an d with f = 1a ■ Noticing that 



J Pt{l(o,ac))vid^ = 'Pu 1 {To>t){vi,l) 



and J PtlArjidfi = ¥ Vl (Xt £ A, Tq > t) (771 , we have shown that v\ is a q.s.d.. 
The rest of the proof is divided into two cases. First assume that /j is a bounded measure. 
Thanks to Theorem l2.3l we know that for any x > 0, any set A C (0, 00) such that 1a £ ¥?(fi) 
and for any t > 1 

F x (X t eA,T >t) = J F y (X t -i e A, T > t — 1) r(l, x, y) n(dy) 

P t -i(lA)(y)r(l,x,y)fi(dy) 

l A (y)(Pt-ir(l,x,.))(y)^dy). 
Since both 1a and r(l, x, .) are in L 2 (/x) and since (H2) is satisfied, we obtain using Theorem 
lira e x ^ t -VF x (X t eA,T >t) = (lA,Vi)n(r(l,x,.),r ll )^. (5.5) 

Since 

-A 



r(l,x,y)rji(y)fj,(dy) = (Pirjx)(x) = e 1 771(0;) 
we get that v\ is the Yaglom limit. 

If fj, is not bounded (i.e. l(o,oo) ^ L 2 (/j)) we need an additional result to obtain the Yaglom 
limit. 

Lemma 5.3. Assume rji £ L 1 (/j) then for all x > 0, there exists a locally bounded function 
Q(x) such that for all y > and all t > 1, 

r(t,x,y) < G(x)e~ Al *r/i(y) . (5.6) 

We postpone the proof of the lemma and indicate how it is used to conclude the proof of the 
theorem. 

If (|5.6p holds, for t > 1, e Xlt r(t,x, .) £ L 1 (/u) and is dominated by ®(x)rji. Since r(l,x, •) £ 
L 2 (/_i) by Theorem E31 using Theorem 13.21 and writing again r(t,x, .) = Pt-\r(l,x, .) fj, a.s., 
we deduce that lim^oo e Xlt r(t,x, .) exists in L 2 (/u) and is equal to 

e Al (r(l,x, .) 5 77i) M r/i(.) = m(x)r]i(-) ■ 

Recall that convergence in L 2 implies almost sure convergence along subsequences. Therefore, 
for any sequence t n — > 00 there exists a subsequence t' n such that 

lim e Xlt ™ r(t' n , x, y) = rji (x) 771 (y) for /^-almost all y > . 
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Since 

F X (T > t' n ) = I r(t' n ,x,y)fj,(dy) , 
Jo 

Lebesgue bounded convergence theorem yields 

/•oo 

lim e Al '- F X (T > t' n ) = m (x) / T ? i(y)/i(dy) . 

That is (|5.5p holds with A = (0, oo) for the sequence t' n . Since the limit does not depend 
on the subsequence, lim^oo e Xlt F x (To > t) exists and is equal to the previous limit, hence 
(|5.5p is still true. The rest of this part follows as before. 

For the last part of the theorem, that is passing from the initial Dirac measures at every 
fixed x > to the compactly supported case, we just use that 6(«) is bounded on compact 
sets included in (0,oo). 

Proof of Lemma \5.3[ According to the parabolic Harnack's inequality (see for example |36j), 
for all x > 0, one can find Qo(x) > 0, which is locally bounded, such that for all t > 1, y > 
and z with \z — x\ < p(x) = h A x 



It follows that 

r(t,x,y) 



2 ' 1 4 

r(t,x,y) < Q (x)r(t + l,z,y). 
{!\z~x\< P (x) r i^ x ^)^( z )^ d < 



\z-x\<p(x) m&Mdz 
, s (J\z-x\<p(x) r ( t + 1 > z iy) 7 ii{z)Kdz) 

< @0{X) 7 

W*-x|<p(x) Vii^Kdz 
(Jr{t+ l,z,y)rji(z)n(dz)) 



< e (x) 

< e (x) 



(l\z-x\<p(x) 

g-Al(t+l) ^y-j 

J\z~x\<p(x) m(z)Kdz] 



since Pt+\f]\ = e x n t + 1 >r] 1 . But Gi(x) = ji^UpM r /i(- z )i u (^ z ) > 0, otherwise r/i, which is a 

solution of the linear o.d.e. \g" — qg' + Aig = on (0, oo), would vanish on the whole interval 
\z — x\ < p(x), hence on (0, oo) according to the uniqueness theorem for linear o.d.e's. The 
proof of the lemma is thus completed with = e _Al Go/@i- □ 



The positive real number Ai is the natural killing rate of the process. Indeed, the limit (15, ip 
obtained in Theorem 15.21 shows that for any x > and any t > 0, 

F x (T >t + 8) = e _ Xlt 
s^oo F X (T > s) 

Let us also remark that 



(T >t) = e- Al *. 
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In order to control the speed of convergence to the Yaglom limit, we first establish the 
following lemma. 

Lemma 5.4. Under conditions (H2) and (H4), the operator P\ is bounded from L°°( / u) to 
L 2 (/x). Moreover, for any compact subset K of (0, oo), there is a constant Ck such that for 
any function f G L 1 ( / u) with support in K we have 



l-Pi/Hi^) < Ck ||/||li( m ) 



Proof. Let g G L°°(/x), since 



we get from (H4) 



\Pig\ < Pi\g\ < lbllL°°( M ) 



\Pig\ dfi < \\g n2 



l°°(m) 



We now recall that (see Section 3) 

p ig (x) = e Q ^ 2 P! (e- Q/2 g\ (x) . 
It follows from Lemma 14.51 that uniformly in x G (0, 1] we have (using hypothesis (H4)) 



Pi [e- Q/2 g) (x) < 0(1) J e-QW* e -i?/iy dy < 

This implies 

J |P l5 |% = J \P X (e-Wg) (x)\ 2 dx < 0(l)\\g\\l^ M , 

and the first part of the lemma follows. For the second part, we have from the Gaussian 
bound of Lemma 14.51 that for any x > and for any / integrable and with support in K 



Pi(e- Q/2 f) (x)\<0(l)J , 



e - Q{y)/2 e-^ 2 / 2 \f(y)\dy 



< 0(1) sup e Q(z)/2 sup e 

z&K z&K 



-(x-z) 2 /2 



~ Qiy) \f(y)\dy<<D(l) e 



x 2 /4, 



K 



- Q{y) \f{y)\dy 



K 



since K is compact. This implies 



J |Pi/|% = J \P X [e- Q ' 2 f) {x)\ 2 dx < Oil) 



2 

Li( M ) • 



□ 



We can now use the spectral decomposition of r(l,x, .) to obtain the following convergence 
result. 

Proposition 5.5. Under conditions (H2) and (H4), for all x > and any measurable subset 
A of (0, oo), we have 



lim e ^- x ^ (F x {X t eA\T > t)-^)) 



(l,Vi) 2 
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Proof. Let h be a non negative bounded function, with compact support in (0, oo). By using 
the semigroup property, Lemma 15,41 and the spectral decomposition for compact self adjoint 
semigroups (see Theorem 13. 2 1) . we have for any t > 2, 

J F x (X t £A,T > t)h{x)dx = (he Q , P t l A )„ = (Pi(/ie Q ) , P^PiU)^ 

= (Pi{heQ) , ( m , PilA) M e- Al (*- 2 ) + (Pi{he Q ) , m ) „ , Pal^e"^'- 2 ) + R(h,A,t) 
with 



j>3 



lL2(/x) 



due to Ai < A2 < A3 < . . ., the Cauchy-Schwarz inequality and Parseval's identity. Note 
that since Pi is symmetric with respect to the scalar product, we have (Pi(he®) , 771}^ = 
e~ Al (he®, rj\)^ and similarly for 772. We also have (771 , P\1a)^ = e~ Al (??i , 1a) /1 and similarly 
for 772. It follows immediately from Lemma [5.4l that for any fixed compact subset K of (0, 00), 
any A and any h satisfying the hypothesis of the proposition with support contained in K, 

\R(h,A,t)\ < 0(l)e- x ^\\he^\\ LHfl) < 0{l)e- x ^\\h\\ U[dx) , 

since h has compact support in (0,oo). Therefore, letting h tend to a Dirac mass, we obtain 
that for any compact subset K of (0, 00), there is a constant Dk such that for any x £ K, 
for any measurable subset A of (0, 00), and for any t > 2, we have 



(X t G A, T > t) - e^ m (x) ( m , Ia)^-^ - e^rj 2 (x) ( m , Ia)^ 1 
The proposition follows at once from 

F x {X t e A,T >t) 



< D K e 



-A 3 t 



,(X t £ A\T >t) 



\{X t e (0,00), T >t) ' 

□ 



6. The Q-process 

As in [5 J (Theorem B), we can also describe the law of the process conditioned to be never 
extinct, usually called the Q-process (also see [2"2"]). 



Corollary 6.1. Assume (H) holds. For all x > and s > we have 

lim ¥ X (X £ B \ Tq > t) = Q X (B) for all B Borel measurable subsets of C([0,s]), 

where <Q X is the law of a diffusion process on (0, 00), with transition probability densities 
(w.r.t. Lebesgue measure) given by 

q (s,x,y)= e ^ s ^-r(s,x,y)e-^ , 
771 (x) 

that is, Q x is locally absolutely continuous w.r.t. ¥ x and 



,{X GB) = E I 1 B {X) e MS -^Yr > r o > s 

771 (x) 
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Proof. First check thanks to Fubini's theorem and k(0 + ) < oo in Hypothesis (HI), that 
A(0 + ) > — oo. We can thus slightly change the notation (for this proof only) and define A as 
A(x) = Jq e®( y } dy. From standard diffusion theory, (A(X tA T );t > 0) is a local martingale, 
from which it is easy to derive that for any y > x > 0, W y (T x < Tq) = A(y)/A(x). 

Now define v(t,x) = p^p-°>*] ■ As in [5, proof of Theorem B], one can prove for any x > 1, 
using the strong Markov property at T x of the diffusion X starting from 1, that v(t,x) < 
A(x)/A(l). On the other hand for x < 1 we obtain v(t,x) < 1. Then for any x > we get 
v(jt,x) < 1 + A(x)/A(l). 

Now thanks to TheoremE21 for all x, e Xlt F x (T > t) -> r/i(x)(l, rji)fj, as t — > oo, and 



hm v[t,x) = — -— . 

t-oo 7/i (1) 

Using the Markov property, it is easily seen that for large t, 

P X (X eB\T >t)=E x [l B (X)v(t - s,X s ),T > s] 



i(T > t 



x(T > t) 



The random variable in the expectation is (positive and) bounded from above by 1 + 
A(X S )/A(1), which is integrable (see below), so we obtain the desired result using Lebesgue 
bounded convergence theorem. 

To see that E x (A(X S ) 1s<t ) is finite, it is enough to use Ito's formula with the harmonic 
function A up to time TqATm- Since A is nonnegative it easily yields K x (A(X S ) 1 s <t at m ) < 
A(x) for all M > 0. Letting M go to infinity the indicator converges almost surely to 1<j<t 
(thanks to Hypothesis (HI)) so the monotone convergence theorem yields (A(X S ) 1 s <t ) < 
A(x). □ 



Recall that v\ is the Yaglom limit. 

Corollary 6.2. Assume (H) holds. Then for any Borel subset B C (0, oo) and any x, 

lim Q X (X S £B)= [ r,l{y)ii{dy) =(??i,l) M f m(y>i{dy). 
s ^°° Jb Jb 

Proof. We know from the proof of Theorem 15.21 that e AlS r(s,x, .) converges to 7?i(x)r?i(.) in 
L 2 (//) as s — ► oo. Hence, since IbVi G ^(a 4 ); 

m (x)Q x (X s G B) = J l B (y)vi{y)e XlS r(s,x,y)fi(dy) m (x) J rfi(y)ii(dy) 

as s — > oo. We remind the reader that dv\ = rjid/j,/ (rji, □ 

Remark 6.3. The stationary measure of the Q-process is absolutely continuous w.r.t. fx, 
with Radon-Nikodym derivative (rji, l)u,Vi which thanks to Proposition 14. II is nondecreasing. 
In particular, the ergodic measure of the Q-process dominates stochastically the Yaglom 
limit. We refer to \27\ [22] for further discussion of the relationship between q.s.d. and 
ergodic measure of the Q-process. 
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7. Infinity is an entrance boundary and uniqueness of q.s.d. 
We start with the notion of quasi limiting distribution q.l.d.. 

Definition 7.1. A probability measure tt supported on (0, oo) is a q.l.d. if there exists a 
probability measure v such that the following limit exists in distribution 

lim ¥ u {X t G • | T > t) = tt(«). 

t— >oo 

We also say that v is attracted to tt, or is in the domain of attraction of it, for the conditional 
evolution. 

Obviously every q.s.d. is a q.l.d., because such measures are fixed points for the conditional 
evolution. We prove that the reciprocal is also true, so both concepts coincide. 

Lemma 7.2. Let tt a probability measure supported on (0, oo). If tt is a q.l.d. then tt is a 
q.s.d.. In particular there exists a > such that for all s > 

P W (T >s) = e 



-as 



Proof. By hypothesis there exists a probability measure v such that lim^oo ¥ u (Xt € • | To > 
t) = 7r(»), in distribution. That is, for all continuous and bounded functions / we have 

} 1U1 TT]) (rp Z~ J.\ = / f(x)ndx). 

t^oo P U (T > t) J 

If we take f{x) = F X (X S G A, Tq > s), since f(x) = f A r(t,x,y)fj,(dy), an application of 
Harnack's inequality and of the dominated convergence theorem ensures that / is continuous 
in (0, oo). 

First, take A = (0, oo), so that f{x) = P x (To > s). Then, we obtain for all s > 

} im p it ^ T\ = F ^ T ° > s >- 
t-^oo r u {iQ > t) 

The left hand side is easily seen to be exponential in s and then there exists a > such that 

¥ n (T >s) = e~ as . 
Second, take f{x) = f x (X s G A,Tq > s) to conclude that 
F n (X s eA,T >s) = lim F u (f(X t )\T > t) = lim F v (X t+s £ A\T >t + s ffB^f 
= e- as 7r(A), 

and then ir is a q.s.d.. □ 

Recall from Theorem 15.21 that under Hypothesis (H), the measure dv\ = rjidfi/(rji, 1) M is the 
Yaglom limit, which in addition is a q.l.d. attracting all initial distribution with compact 
support on (0, oo). It is natural to ask about the uniqueness of the q.s.d.. Here again, 
our assumptions on the behavior of q at infinity will allow us to characterize the domain of 
attraction of the q.s.d. v\ associated to r\\. This turns out to be entirely different from the 
cases studied in [5] for instance. 

We say that the diffusion process X comes down from infinity if there is y > and a time 
t > such that 

limPa,(r B < t) > 0. 

xfoo 
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This terminology is equivalent to the property that oo is an entrance boundary for X (for 
instance see [31] page 283]). 

Let us introduce the following condition 

/OO /'OO 
e Q(y) J e ~ Q{z) dzdy < oo. 

Tonelli's Theorem ensures that (H5) is equivalent to 

/oo py 
e -Q(y) J e Q(*)dz dy < oo. (7.1) 

If (H5) holds then for y > 1, J°° e~ Q ^ dz < oo. Applying the Cauchy-Schwarz inequality 

we get (x — l) 2 = (J** e Q / 2 e~ Q / 2 dz) 2 < J® e®dz e~ Q dz, and therefore (H5) implies that 
A(oo) = oo. 

Now we state the main result of this section. 

Theorem 7.3. Assume (H) holds. Then the following are equivalent: 

(i) X comes down from infinity; 

(ii) (H5) holds; 

(hi) v\ attracts all initial distributions v supported in (0, oo), that is 

lim F v (X t e • | T > t) = 

In particular any of these three conditions implies that there is a unique q.s.d.. 

Remark 7.4. It is not obvious when Condition (H5) holds. In this direction, the following 
explicit conditions on q, all together, are sufficient for (H5) to hold: 

• q{x) > qo > for all x > xq 

• limsup x ^ O0 q'{x)/2q 2 (x) < 1 

• r°° -rr dx < oo . 

Jx q(x) 

Indeed, check first that these conditions imply that q{x) goes to infinity as x — * oo. Then 
defining s(y) := e~®^ dz, the first condition above implies that s(y)e®^ is bounded in 
y > xq. Integrating by parts on J se^dz gives 

rx rx S x f' X 1 f X a' 

/ se Q dz = / — e Q 2qdz = — e Q + / — dz + / se Q -^dz. 
Jxq •> XO 2q 2g Ixo J XQ 2q J XQ 2q z 

Since se® /2q vanishes at infinity, the third condition implies that se^{l — q'/2q 2 ) is integrable 
and thanks to the second condition we conclude that (H5) holds. 

On the other hand, if (H5) holds, q'{x) > for x > xq and q(xo) > 0, then q(x) goes to 
infinity as x — > oo and dx < oo . 

We can retain that under the assumption that q'(x) > for x > xq and q(x) goes to infinity 
as x — > oo, then 

f°° 1 

(H5) / — — dx < oo . 

Ji q{ x ) 



24 P. CATTIAUX, P. COLLET, A. LAMBERT, S. MARTINEZ, S. MELEARD, AND J. SAN MARTIN 

Indeed, the only thing left to prove is the sufficiency of (H5). Since s(y) tends to as y —* oo 
(because Q grows at least linearly), then by the mean value theorem we have 

f°° e-«Mcfe l 
Jy 1 



where £ G [y, oo). Using that g is monotone we obtain the bound 

f°° e-«Wd2 i 



e -Q(v) ~ 2q{yY 

and the equivalence is shown. <0> 
The proof of Theorem 17.31 follows from Propositions 17.51 17.61 and 17.71 

Proposition 7.5. Assume (HI) holds. If there is a unique q.s.d. that attracts all initial 
distributions supported in (0, oo), then X comes down from infinity. 

Proof. Let ir be the unique q.s.d. that attracts all distributions. We know that W n (To > t) = 
e —at £ or some q, > o. Since absorption is certain then a > 0. For the rest of the proof let 
v be any initial distribution supported on (0, oo), which by hypothesis is in the domain of 
attraction of ir that is, for any bounded and continuous function / we have 

/■oo roo 

lim / F u (X t £dx\T > t)f{x) = / /(x)Tr(dx). 
Jo Jo 

We now prove that for any A < a, E^(e AT °) < oo. As in Lemma 17.21 we have for any s 

i im p - (r ° >t+ '> =e -°-. 

t-oo P„(T > t) 

Now pick A G (0, a) and e > such that (1 + e)e A_Q < 1. An elementary induction shows 
that there is to such that for any t > to, and any integer n 

(To > t + 



IP, (To > t) 

Breaking down the integral J* i ^°P l/ (To > s) e Xs ds over intervals of the form (n, n + 1] and 
using the previous inequality, it is easily seen that this integral converges. This proves that 
^i/(e XT °) < oo for any initial distribution v. 

Now fix A = a/2 and for any x > 0, let g(x) = K x (e XT °) < oo. We want to show that g is 
bounded, which trivially entails that X comes down from infinity. Thanks to the previous 
step, for any nonnegative random variable Y with law v 

E(g(Y)) =E„(e ATo ) < oo. 

Since Y can be any random variable, this implies that g is bounded. Indeed, observe that g 
is increasing and g(0) = 1, so that a := l/g(oo) is well defined in [0, 1). Then check that 

(A \ g '^> A 

vm = n VT\2 dx 

(1 - a)g{x) 2 
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is a probability density on (0, oo). To conclude we use the fact that J gdv < oo to get 

9 J (l-a)g(x) X l- a 9K >\* I- a ' 
and then g is bounded. 

□ 

Proposition 7.6. The following are equivalent 

(i) X comes down from infinity; 

(ii) (H5) holds; 

(hi) for any a > i/iere exists y a > suc/i i/iai sup 3;>2/a [e" 7 *" 1 ] < oo . 

Proof. Since (i) is equivalent to oo being an entrance boundary and (ii) is equivalent to (|7.ip 
we must show that "oo is an entrance boundary" and (|7.ip are equivalent. This will follow 
from [17, Theorem 20.12,(w)]. For that purpose consider Yt = A(Xt). Under each one of 
the conditions (i) or (ii) we have A(oo) = oo. It is direct to prove that Y is in natural scale 
on the interval (A(0), oo), that is, for A(0) <a<y<b<oo = A(oo) 

P„Ctf < T?) = 

where is the hitting time of a for the diffusion Y. Then, oo is an entrance boundary for 
Y if and only if 

/•oo 

/ y m(dy) < oo, 
J o 

where m is the speed measure of Y, which is given by 

m(dy) = (A'(A-W ' 
see [HI formula (5.51)], because Y satisfies the SDE 

dY t = M(K- x (Y t ))dB t . 
After a change of variables we obtain 

i f'OQ FX 

ym(dy)= \ e~ Q{y) / e Q{z) dzdx. 



Therefore we have shown the equivalence between (i) and (ii). 

We continue the proof with (ii) => (Hi). Let a > 0, and pick x a large enough so that 

/•oo poo 1 

J Xa L - 2a 

Let J be the nonnegative increasing function defined on [x ai oo) by 

f'X f'OQ 

J( x ) = / e Q(y) / e - Q{z) dzd y . 



Then check that J" = 2qJ' — 1, so that LJ = —1/2. Set now y a = 1 + x a , and consider a 
large M > x. Ito's formula gives 

Ex(e a(MTMATto) ^(^AT M AT a J) = J(x) +~E X I J e as (aJ(X s ) + LJ(X S )) dsj . 
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But LJ = -1/2, and J(X S ) < J(oo) < l/(2a) for any s <T ya , so that 

Ex[e a(tAT M AT m ) J iXtATMAT J] < J(x) . 

But J is increasing, hence for x > y a one gets l/(2a) > J(x) > J(y a ) > 0. It follows 
that E x {e a( - tATMAT y^) < l/(2aJ(y a )) and finally E x (e aT y«) < l/(2aJ(y a )), by the monotone 
convergence theorem. So (Hi) holds. 

Finally, it is clear that (Hi) =^ (i). □ 

Proposition 7.7. Assume (H) holds. If there is xq such that sup a; > a . ~E x (e XlTx o ) < oo, then 
v\ attracts all initial distribution supported in (0,oo). 

The proof of this result requires the following control near and oo. 

Lemma 7.8. Assume (H) holds, and sup a;>3 , E x (e XlTx o) < oo. For h £ L 1 ( / u) strictly 
positive in (0, oo) we have 

lim lim sun S«m^{T >t)^ dx ) _ 

lim lim sud Im H^(T >tMdx) _ , . 

Proof. We start with (|7.2p . Using Harnack's inequality, we have for e < 1 and large £ 
JoMx^CTq >t)n(dx) F 1 (T >t) ghi^jdz) 



J h(x)F x (T > Cr(t - 1, 1, 1) ft h(x) l x(dx) ft n(dy) ' 
then 

/„ e h(x)F x (T > t)fi(dx) Pi (r > t) J' /i(z)^(d2! 

Inn sup -j — — — < lim sup 



t^oo fh{x)¥ x (To>t)n(dx) t^oo Cr(t- 1,1,1) ft h(x)^(dx) ft fx(dy) 

C r?i(l) J x 2 h(x)n(dx) ft n(dy) ' 
and the first assertion of the lemma is proven. 

For the second limit, we set Aq := sup E x (e XlTx o ) < oo. Then for large M > xq, we have 

x>xo 

(To >t)= I P 20 (T > u)P x (r x . ed(f-u)) + P^T^ > t). 

JO 



Using that lim e Ain P :Co (To > u) = r/i(a;o)(r/i, 1)» we obtain that B := supe XlU F Xo (T > 
u) < oo. Then 



P. (T > t) <B £ e^T^ e d(t - «)) + P^T^ > t) 

< 5 e- Al *E a; (e AlTa; o) + e - Ali E x (e AlT:c o) < e~ Ali A)(£o + 1), 

and (|7.3p follows immediately. □ 
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Proof of Proposition \7. 1\ Let v be any fixed probability distribution whose support is con- 
tained in (0,oo). We must show that the conditional evolution of v converges to v\. We 
begin by claiming that v can be assumed to have a strictly positive density h, with respect 
to fi. Indeed, let 



r(l, x, y)v(dx). 



Using Tonelli's theorem we have 
r(l,x,y)v(dx) n(dy) = 



'(l,x,y)iM(dy)u(dx) = J F X (T > l)u(dx) 



< 1, 



which implies that j r(l,x,y)v(dx) is finite dy—a.s.. Also I is strictly positive by Harnack's 
inequality. Finally, define h = £/ j idjjL. Notice that for dp = hd/j, 

M X t+i G • I T > t + 1) = f P (X t e • I T > t), 

showing the claim. 

Consider M > e > and any Borel set A included in (0, 00). Then 

jP x (X t £A,T > t)h(x)n{dx) f e M W x (X t eA,T > t)h(x)n(<k 



fP x {T > t)h(x)n(dx) J e M V x (T > t)h{x)p{dx) 

is bounded by the sum of the following two terms 



II 



12 



fF x (X t £A,To>t)h(x)n(dx) _ J e M P x (X t £A, T >t)h(x) fj,(dx) 
f¥ x {T >t)h{x)n{dx) f¥ x {T >t)h(x) n(dx) 

J g M f x {X t dA, T >t)h(x) fj,(dx) J c M P x (X t £A, T >t)h(x) fj,(dx) 



f¥ x {T >t)h(x)fi{dx) 



C P x (T >t)h(x) fi(dx) 



We have the bound 



II V 12 < 



/ o e P«(T > t)h(x)fi(dx) + J^Pz(T > t)h(x) n{dx) 



JP x (T >t)h(x)p{dx) 
Thus, from Lemma 17.81 we get 

jP x (X t £ A, T > t)h(x)(i{dx) f € M Px(X t eA,T > t)h(x) fi{dx) 



lim lim sup 

40,MToo t -+oc 



JF X (T > t)h(x)fi(dx) 



jf ¥ X (T > t)h(x) n(dx) 



On the other hand we have using 



J e M P x (X t e A, T > t)h(x)p(dx) _ f A rn(z)n(dz) 



lim 

t— >oo 



f e M F X (T > t)h(x) n(dx) " Jr+ m{z)n(dz) 
independently of M > e > 0, and the result follows. 



□ 



The following corollary of Proposition 17.61 describes how fast the process comes down from 
infinity. 

Corollary 7.9. Assume (H) and (H5) hold. Then for all A < Ai, sup a . >0 E x [e AT °] < 00. 
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Proof. We have seen in Section [51 (Theorem I5.2p that for all x > 0, lim^oo e Xlt F x (To > t) = 
r]i(x)(r]i,l) fl < oo i.e. E x [e ATo ] < oo for all A < Ai. Applying Proposition 17.61 with a = A 
and the strong Markov property it follows that sup x>yx M x [e ) ' To ] < oo. Furthermore, thanks 
to the uniqueness of the solution of (jl.ip . Xf < X\ x a.s. for all t > and all x < y\, hence 
E x [e A:ro ] < Ey A [e AT °] for those x, completing the proof. □ 

The previous corollary states that the killing time for the process starting from infinity has 
exponential moments up to order Ai. In [21] an explicit calculation of the law of To is done 
in the case of the logistic Feller diffusion Z (hence the corresponding X) and also for other 
related models. In particular it is shown in Corollary 3.10 therein, that the absorption time 
for the process starting from infinity has a finite expectation. As we remarked in studying 
examples, a very general family of diffusion processes (including the logistic one) satisfy all 
assumptions in Corollary 17.91 which is thus an improvement of the quoted result. 

We end this section by gathering some known results on birth-death processes that are close 
to our findings. Let Y be a birth-death process with birth rate A n and death rate fj, n when 
in state n. Assume that Ao = no = and that extinction (absorption at 0) occurs with 
probability 1. Let 

S = TTi + ^(A n 7T n ) _1 TTi , 

i>l n>l i>n+l 

where 

A1A2 • • • A n _! 

vr n = • 

fJ-lfJ-2 ■■■ 

In this context sure absorption at 0, that is (HI), is equivalent to A := ^j> 1 (Aj7Tj) _1 = 00 
(see pjJJ formula 7.9]). On the other hand we also have Ei(To) = Yj%>i ^i- We may state 

Proposition 7.10. For a birth-death process Y that satisfies (HI), the following are equiv- 
alent: 

(i) Y comes down from infinity; 

(ii) There is one and only one q.s.d.; 

(iii) lim nToo j E n (T ) < 00; 

(iv) S < 00. 

Proof. In |8j Theorem 3.2], it is stated the following key alternative: S < 00 iff there is a 
unique q.s.d.; if S = 00 then there is no q.s.d. or there are infinitely many ones. Then the 
equivalence between (ii) and (iv) is immediate. Also it is well known that 

n-l 

E„(T ) = ^7Tj + ^(A^)" 1 ^ ^i, 

i>l r=l «>r+l 

see for example |19[ formula 7.10]. Therefore (iii) and (iv) are equivalent. Let us now examine 
how this criterion is related to the nature of the boundary at +00. From the table in [TJ 
section 8.1] we have that +00 is an entrance boundary iff A = 00, Ei(To) < 00 and S < 00. 
Finally, S < 00 implies that Ei(To) < 00 and this ensures Pi (To < 00) = 1, that is A = 00. 
This shows the result. □ 
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8. Biological models 

8.1. Population dynamics and quasi-stationary distributions. Our aim is to model 
the dynamics of an isolated population by a diffusion Z := {Zt\t > 0). Since competition 
for limited resources impedes natural populations with no immigration to grow indefinitely, 
they are all doomed to become extinct at some finite time To. However, To can be large 
compared to human timescale and it is common that population sizes fluctuate for large 
amount of time before extinction actually occurs. The notion of quasi-stationarity captures 
this behavior [301134] , 

The diffusions we consider arise as scaling limits of general birth-death processes. More 
precisely, let (Z n )n^ be a sequence of continuous time birth-death processes Z N := 
{Z^;t > 0), renormalized by the weight iV -1 , hence taking values in iV _1 N. Assume that 
their birth and death rates from state x are equal to &tv(x) and disr(x), respectively, and 
&jv(0) = ^jv(O) = 0, ensuring that the state is absorbing. We also assume that for each N 
and for some constant Bn, 

b N (x) < (x + 1)B N , x > 
and that there exist a nonnegative constant 7 and a function h £ C 1 ([0, 00)), h(0) = 0, called 
the growth function, such that 

VxG(0,oo): lim -^(b N (x) - d N (x)) = h(x) ; lim 7^p?(b N (x) + d N (x)) = jx. (8.1) 

Important ecological examples include 

(i) The pure branching case, where the individuals give birth and die independently, so 
that one can take &at(x) = (7JV + X)Nx and d/v(x) = (jN+fj,)Nx. Writing r := A — /i 
for the Malthusian growth parameter of the population, one gets h(z) = rz. 

(ii) The logistic branching case, where b^{x) = ('yN + X)Nx and dpf(x) = (jN + fi)Nx + 
jjNx(Nx-l). The quadratic term in the death rate describes the interaction between 
individuals. The number of individuals is of order N, the biomass of each individual 
is of order N , and c/N is the interaction coefficient. The growth function is then 
h(z) = rz — cz 2 . 

(iii) Dynamics featuring Allee effect, that is, a positive density-dependence for certain 
ranges of density, corresponding to cooperation in natural populations. A classical 
type of growth function in that setting is h(z) = rz(j^ — 1)(1 — Observe that in 
this last case, the individual growth rate is no longer a monotone decreasing function 
of the population size. 

Assuming further that (Zq^tvsn converges as ./V — > 00 (we thus model the dynamics of a 
population whose size is of order N), we may prove, following Lipow [24] or using the tech- 
niques of Joffe-Metivier [16], that the sequence (Z n )n^ converges weakly to a continuous 
limit Z. The parameter 7 can be interpreted as a demographic parameter describing the 
ecological timescale. There is a main qualitative difference depending on whether 7 = or 
not. 

If 7 = 0, then the limit Z is a deterministic solution to the dynamical system Zf = h(Zt). 
Since h(0) = 0, the state is always an equilibrium, but it can be unstable. For example, 
in the logistic case h(z) = rz — cz 2 and it is easily checked that when r > 0, the previous 
dynamical system has two equilibria, which is unstable, and r/c (called carrying capacity) 
which is asymptotically stable. In the Allee effect case, and K are both stable equilibria, 
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but Kq is an unstable equilibrium, which means the population size has a threshold Kq to 
growth, below which it cannot take over. 

If 7 > 0, the sequence (Z n )n^ converges in law to the process Z, solution to the following 
stochastic differential equation 

dZ t = y/iZtdBt + h{Z t )dt. (8.2) 

The acceleration of the ecological process has generated the white noise. Note that /i'(0 + ) is 
the mean per capita growth rate for small populations. The fact that it is finite is mathe- 
matically convenient, and biologically reasonable. Since h(0) = 0, the population undergoes 
no immigration, so that is an absorbing state. One can easily check that when time goes 
to infinity, either Z goes to oo or is absorbed at 0. 

When h = 0, we get the classical Feller diffusion, so we call generalized Feller diffusions the 
diffusions driven by (|8.2p . When h is linear, we get the general continuous-state branching 
process with continuous paths, sometimes also called Feller diffusion by extension. When h 
is concave quadratic, we get the logistic Feller diffusion [9| \2~T\. 

Definition 8.1. (HH) We say that h satisfies the condition (HH) if 

r\ v h ( x ) r -\ v xh '( x ) n 

(i) lim — — = — oo, (n) lim = U. 

In particular (HH) holds for any subcritical branching diffusion, and any logistic Feller diffu- 
sion. Concerning assumption (i), the fact that h goes to — oo indicates strong competition in 
large populations resulting in negative growth rates (as in the logistic case) . Assumption (ii) 
is fulfilled for most classical biological models, and it appears as a mere technical condition. 

Gathering all results of the present paper and applying them to our biological model yields 
the following statement, which will be proved at the end of this section. 

Theorem 8.2. Let Z be the solution of 118. We assume h £ C 1 ([0, oo)), h(0) = and 
that h satisfies assumption (HH). Then, for all initial laws with bounded support, the law of 
Zt conditioned on {Zf ^ 0} converges exponentially fast to a probability measure v , called the 
Yaglom limit. The law Q_ x of the process Z starting from x and conditioned to be never extinct 
exists and defines the so-called Q-process. This process converges, as t — > oo, in distribution, 
to its unique invariant probability measure. This probability measure is absolutely continuous 
w.r.t. v with a nondecreasing Radon-Nikodym derivative. 
In addition, if the following integrability condition is satisfied 

dx 

^) <0 °' 

then Z comes down from infinity and the convergence of the conditional one- dimensional 
distributions holds for all initial laws. In particular, the Yaglom limit v is then the unique 
quasi- stationary distribution. 



Proof. For Z solution to (pT2j) . recall that X t = 2y / Z t /j satisfies the SDE dX t = dB t -q(X t )dt 
with 

. . 1 2h(jx 2 /A) 
q( x ) = ii — <—L x> 0. 

2x 7X 
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In particular we have q'(x) = —-^p + 2h ^ c J / ^ > — h'(jx 2 /A) and 

q 2 (x) - q'(x) = A + (A^( 7 * 2 /4) - A) + ^V*)- 

Under assumption (HH) we have the following behaviors at and oo: q(x) ~ l/2x, as well 

xj.0 

as 

^2 



.r 



i{x) ~ A and ( , 2 _ ,o(2v^) ~ MA (I + A^) 



Then, it is direct to check that hypothesis (H2) holds 



lim q 2 (x) — q'(x) = oo , C := — inf q 2 (x) — q'(x) < oo. 

x-^oo xe(0,oo) 



We recall that 

Q(x) = 2q{y)dy, A(z) = £ e Q ^dy and k(x) = e Q ^ ^V Q ^dzj dy. 

Straightforward calculations show that 

lim = oo and A := lim (Q(x) — log(x)) G (—00,00). 

In particular, A(oo) = 00 and the integrand in the definition of k is equivalent to ylog(y) 
which ensures k(0 + ) < 00. Thus X, and consequently Z, is absorbed at with probability 
1, that is hypothesis (HI) holds. 
We now continue with (H3) which is 

' e- Q{y) dy < 00. 



q 2 {y)-q'{y) + C + 2 
This hypothesis holds because near the integrand is of the order 

— e W(y > — y. 

W 

For the first part of the Theorem it remains only to show that (H4) holds 

00 pi 



e W( >dx < 00 and / x e ^ { " dx < 00. 
'1 Jo 

The first integral is finite because Q grows at least linearly near 00 and the second one is 

finite because the integrand is of order l/y/x for x near 0. 

Hence we can apply Theorem 15.21 Proposition 15.51 an d Corollaries 16.11 and 16.21 to finish with 
the proof of the first part of the Theorem. 

For the last part of the Theorem we need to show that X comes down from infinity which 
is equivalent to (H5) . Thanks to Remark 17.41 there is a simple sufficient condition for this 
hypothesis to hold, which has three components. The first one 

q{x) > go > for all x> xq 

follows from (HH)(z). The second one 

limsupg / (x)/2g 2 (x) < 1 
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is equivalent to 



xh'(x) 2 



which clearly follows from (HH)(m). Finally the third one 

POO ^ 



■ dx < oo, 



thanks to (HH)(i), is equivalent to 



dx = [ — ; . , dz < oo. 



2/i( 7 x 2 /4) J -h(z) 
This is exactly the extra assumption made in the Theorem and the result is proven. □ 

8.2. The growth function and conditioning. Referring to the previous construction 
of the generalized Feller diffusion (|8.2j) . we saw why h(z) could be viewed as the expected 
growth rate of a population of size z and h{z) / z as the mean per capita growth rate. Indeed, 
h(z) informs of the resulting action of density upon the growth of the population, and h(z)/z 
indicates the resulting action of density upon each individual. In the range of densities z 
where h(z)/z increases with z, the most important interactions are of the cooperative type, 
one speaks of positive density- dependence. On the contrary, when h(z)/z decreases with z, 
the interactions are of the competitive type, and density-dependence is said to be negative. 
In many cases, such as the logistic one, the limitation of resources forces harsh competition 
in large populations, so that, as z — ► oo, h(z)/z is negative and decreasing. In particular h(z) 
goes to — oo. The shape of h at infinity determines the long time behavior of the diffusion Z. 

Actually, if h goes to infinity at infinity, such as in the pure branching process case (where 
h is linear), Theorem 18.21 still holds if (HH)(i) is replaced with the more general condi- 
tion lim^oo = dboo, provided the generalized Feller diffusion is further conditioned on 
eventual extinction. Indeed, the following statement ensures that conditioning on extinction 
roughly amounts to replacing h with —h. 

Proposition 8.3. assume that Z is given by \8.2i) . where h £ C 1 ([0, oo)), h(0) = 0, 
linx^oo = oo. Define u(x) := P 2 ( lim Z% = 0) and let Y be the diffusion Z conditioned 
on eventual extinction. Then Y is the solution of the SDE, Yq = Zq 

dY t = ^Y t dB t + (h(Y t ) + jY t ^l) dt. (8.3) 
If, in addition h satisfies (HH)(ii) then 

1 1 \ u '(y) 1 1 \ 

Kv) + jy—rr ~2/^oo -h{y). 
u{y) 

Proof. Let J(x) := Jq dz which is well-defined since h G C 1 ([0, oo)) with h(0) = 0. We 
set 

roo 

v(x) := a / e~ J{z) dz, 
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with a = (J °° e~ J ( z 'dz)~ 1 (well-defined by the growth of h near oo). Now we prove that 
u = v. It is easily checked that v is decreasing with v(0) = 1, i>(oo) = 0, and that it satisfies 
the equation ^xv"(x) + h(x)v'(x) = for all x > 0. 

As a consequence, (v(Z t );t > 0) is a (bounded hence) uniformly integrable martingale, so 
that 

v{x) = E x (v(Z t )) -> v(0)F x ( lim Z t = 0) + v(oo)F x ( lim Z t = oo) = u(x), 
as i — > oo, so that indeed u = v. 

Using the strong Markov Property of Z we obtain that for any Borel set A C (0, oo) and 
s > 

P.(y. € A) = P.(Z. € A | To < oo) — E, ( p*ff * g . Z. € a) = Ej (^, Z. 6 a) ■ 
Then for any measurable and bounded function / we get 

Now if / is C 2 and has compact support contained in (0, oo), we get from Ito's formula that 
uf is in the domain of L z , the generator of and then / is in the domain of the generator 
L Y of Y and moreover 

L Y (f)(x) = -±-L z (uf)(x) = Zxf(x) + (h{x) +1 x^\ 
u[x) I \ u(x) 

Then, since h is locally Lipschitz we obtain that the law of Y is the unique solution to the 
SDE (ED. 



Let us show the last part of the proposition. Notice that J is strictly increasing after some 
xq, so we consider its inverse (f on [J(xq), oo). Next observe that for x > xq, 



oo ;>oo 



' 7 (x) = e J & e~ J ( z) dz = e J{x) e- b <p'(b)db, 

Jx Jj(x) 



with the change b = J{z). As a consequence, we can write for y > J(xq) 

ii f'°° r°° 

- -(cp(y)) = e y / e-y (6)d& = / e-y (y + b)db. (8.4) 



o 



Because /i tends to oo, J(x) > (1+s) log(x) for x sufficiently large, so that (f(y) < exp(y/(l + 
e)), and <£>(y)exp(— y) vanishes as y — > oo. Now, since 

<p'(y) exp(— y) also vanishes. Since h is differ entiable, J is twice differentiable, and so is tp, 
so performing an integration by parts yields 

<p'(y)= e- b <p'{y + b)db - e- b ip"(y + b)db. (8.5) 
Jo Jo 

Since <//(J(x)) = 1/ J'(x), we get J'(x)ip"(J(x)) = (1/J'(x))', so by the technical assumption 

(HH)(«), 
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as x — ► oo. Then, as y — > oo we have ip"(y) = o(<p'(y)). This shows, thanks to (|8,5p . that 

/•oo 

/ e-V(l/ + 6)*~»-oo <f/(y) 
Jo 

which entails, thanks to (|8.4|) . that 

J My)) ~j/^oo ¥>'(y). 

This is equivalent to 

u' 

7x— (x) ~ x ->oc -7xJ'(x) = -2h(x), 
u 

which ends the proof. □ 

Let us examine the case of the Feller diffusion (pure branching process), where h(z) = rz. 
First, it is known (see e.g. [231 Chapter 2]) that when r > 0, the supercritical Feller diffusion 
Z conditioned on extinction is exactly the subcritical Feller diffusion with h(z) = —rz. The 
previous statement can thus be seen as an extension of this duality to more general population 
diffusion processes. 

Second, in the critical case (r = 0), our present results do not apply. Actually, the (critical) 
Feller diffusion has no q.s.d. [22]. Third, in the subcritical case (r < 0), our results do apply, 
so there is a Yaglom limit and a Q-process, but in contrast to the case when 1/h is integrable 
at oo, it is shown in [22] that subcritical Feller diffusions have infinitely many q.s.d.. 



Appendix A. Proof of Lemma 14.51 



We first prove the second bound. For any nonnegative and continuous function / with 
support in M + we have from hypothesis (H2) 

(q 2 ~ q'){u s )ds 



pi(x,u) f{u)du = E 



/Ml)) IktqM exp 



1 



<e c / 2 E w M/Ml))li<T H] • 
The estimate (|4.3|) follows by letting f(z)dz tend to the Dirac measure at y with K3 = e c ^ 2 , 
that is 

Pl(x,y) < K 3 pf(x,y). 
Here pf (x,y) = ^= (e~( x ~ y f ' I 2 — e~^ x+y ^ 2 / 2 j (see for example [18, page 97]). 

Let us now prove the upper bound in (|4.2|) . Let B\ be the function defined by 



Bt(z) := inf (q 2 (u) - q'(u)) . 

u>z 



We have 



/ 



Pi(x,y) f(y)dy = E 



/Ml)) 1i<t 1i<t i/3 exp 



{q 2 - q'){u s )ds 



/Ml)) li<To Ii>t« /3 exp \-\ J (q 2 - q')^s)d 
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For the first expectation we have 
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E 



/M 1 )) IktqIkt^ exp 



1 



(q 2 ~ q')(u s )ds 



< e -B!(x/3)/2 gW x [/(^(l)) Il<T(j ] . 

For the second expectation, we obtain 

-l 



/M 1 )) Ii<T *i>t x/3 exp 



1 



(q - q')(u s )ds 



< e C/2 ^ 



f(u(l))l 1<To ll>T x 



/3 



e c/2 r E m x [/Ml)) Ii<Tq] _ E w. r /Ml)) j 



x/3 



Using a limiting argument as above and the invariance by translation of the law of the 
Brownian motion, and firstly assuming that y/2 < x < 2y, we obtain 

Pi(x, y) < e- B ^/ 2 p?(x, y) + e c ' 2 (pf(x, y) - p? (2x/3, y - x/3)) , 



p?(x,y)-p^(2x/3,y-x/3) 



D 



-(y+ x /3) 2 /2 _ e -(x+y)V2\ < 



Since the function B\ is non-decreasing, we get for y/2 < x < 2y 



max{x,y} 2 /18 



Pi{x,y) < 



e -B 1 (ma,x{x,y}/6)/2 , g - max{x,y} 2 /18 



If x/y ^]l/2, 2[, we get from the estimate 



, , K 3 
Pi{x,y) < e 



We now define the function B by 



(y-x) 2 /2 < 



K 3 



■ max{ic,y} 2 /8 



2tt 



2vr 



B( : ) : - Ior [ ) + minlSx^/e)^ , z 2 /36} . 



2tt 

It follows from hypothesis (H2) that lim^oo = oo. Combining the previous estimates 
we get for any x and y in M + 

Pi(x,y) < e - 2B (™*{x>y}) . 

The upper estimate (|4.2p follows by taking the geometric average of this result and A4.3H . We 
now prove that p~\(x,y) > 0. For this purpose, let a = min{x,y}/2 and b = 2max{x,y}. We 
have as above for every nonnegative continuous function / with support in M + 



pi(x,y) f(y)dy>¥} 



/Ml)) I KT [a , i)] exp 

z ./o 



(<7 2 - q'){v s )ds 



where we denote T[ a ,b] the exit time from the interval [a, 6]. Let 

R a ,b = sup (g 2 (x) - q'(x)) , 
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this quantity is finite since q G C 1 ((0, oo)). We obtain immediately 

J Pi(x,y) f(y)dy > e~ R ^ 2 f p [ ? M (x,y) f(y)dy 

where we denote p[ a '^ the heat kernel with Dirichlet conditions in [a, b\. The result follows 
from a limiting argument as above since p^' b \x,y) > 0. □ 
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